x1=ones(1,29570);
x5=ones(1,289)+4;
x2=ones(1,14134)+1;
x3=ones(1,6819)+2;
x4=ones(1,1711)+3;
x6=ones(1,24)+5;
x8=ones(1,3)+7;
x7=ones(1,8)+6;
x9=ones(1,1)+8;
x=[x1,x2,x3,x4,x5,x6,x7,x8,x9];
h1 = histogram(x');
hold on
y10=zeros(1,253)+10;
y13=zeros(1,79)+13;
%y26=zeros(1,10)+26;
y11=zeros(1,178)+11;
y9=zeros(1,355)+9;
y6=zeros(1,1446)+6;
y7=zeros(1,909)+7;
y2=zeros(1,9094)+2;
y5=zeros(1,2300)+5;
%y39=zeros(1,4)+39;
y8=zeros(1,621)+8;
y16=zeros(1,40)+16;
y17=zeros(1,35)+17;
%y28=zeros(1,5)+28;
y4=zeros(1,3515)+4;
%y24=zeros(1,6)+24;
y1=zeros(1,27783)+1;
%y83=zeros(1,1)+83;
y12=zeros(1,122)+12;
y19=zeros(1,17)+19;
%y32=zeros(1,2)+32;
%y40=zeros(1,1)+40;
%y54=zeros(1,2)+54;
%y34=zeros(1,3)+34;
%y30=zeros(1,4)+30;
y3=zeros(1,5554)+3;
y15=zeros(1,40)+15;
%y25=zeros(1,5)+25;
%y52=zeros(1,1)+52;
y14=zeros(1,78)+14;
%y22=zeros(1,14)+22;
%y75=zeros(1,1)+75;
%y21=zeros(1,17)+21;
%y23=zeros(1,11)+23;
%y36=zeros(1,1)+36;
%y41=zeros(1,1)+41;
y18=zeros(1,24)+18;
y20=zeros(1,13)+20;
%y33=zeros(1,2)+33;
%y48=zeros(1,1)+48;
%y27=zeros(1,4)+27;
%y29=zeros(1,4)+29;
%y38=zeros(1,1)+38;
%y37=zeros(1,1)+37;
%y49=zeros(1,1)+49;
y=[y10,y13,y11,y9,y6,y7,y2,y5,y8,y16,y17,y4,y1,y12,y19,y3,y15,y14,y18,y20];
h2 = histogram(y');
xlabel('# of Locations')
ylabel('# of Times')
legend('P\zeta=0.9','P\zeta=1.0')